Bound-Free Electron-Positron Pair Production in Relativistic Heavy-Ion Collisions 
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The bound-free electron-positron pair production is considered for relativistic heavy ion collisions. 
In particular, cross sections are calculated for the pair production with the simultaneous capture of 
the electron into the Is ground state of one of the ions and for energies that are relevant for the 
Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Colliders (LHC). In the framework 
of perturbation theory, we applied Monte-Carlo integration techniques to compute the lowest-order 
Feynman diagrams amplitudes by using Darwin wave functions for the bound states of the elec- 
trons and Sommerfeld-Maue wave functions for the continuum states of the positrons. Calculations 
were performed especially for the collision of Au + Au at 100 GeV/nucleon and Pb + Pb at 3400 
GeV/nucleon. 

PACS numbers: 25.75.Dw; 25.30.Rw 
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I. INTRODUCTION 

The bound-free electron-positron pair production plays 
an important role at modern colliders such as the Rela- 
tivistic Heavy Ion Collider (RHIC) or the Large Hadron 
Colliders (LHC), since it may restrict the luminosity of 
the ion beams that will be available. Especially in pe- 
ripheral collisions of the ions, it is known that a large 
number of lepton pairs can be produced owing to the 
Lorcntz contracted electromagnetic fields that occur in 
course of the collisions. In the bound-free pair produc- 
tion, the electron is captured by one of the colliding ions 



Z a + Zb 



{Z a 



Z b 



(1) 



and leads to the loss of the (one-electron) ion from 
the beam. The bound-free pair production (BFPP) 
will therefore be an important problem at the LHC; in 
fact, this process does not only reduce the intensity of 
the beam but also leads to a separate beam of one- 
electron ions that strikes the beam-pipe about 140 meters 
away from the interaction point. In the worst scenario, 
there might be enough energy in this separated beam to 
quench the LHC magnets as it was first pointed out by 
Bruce et alAM but was investigated in further detail in 
Refs. [II If 

A first computation on the bound-free pair production 
cross sections were performed by Baltz, Rhoadcs-Brown 
and Weneser [f| in the mid 1990ies. These authors used 
large-basis coupled-channel Dirac-equation of BFPP in 
their calculations. In particular, Baltz and coworkers de- 
rived a simple formula for the BFPP cross sections at 
ultra-relativistic energies (7 = 23 000) 



obfpp = Aln (7) + B 



(2) 



where A and B are the parameters independent of en- 
ergy. In this expression for the cross section, the Alnj 
term represents the region of large impact parameters, 
and A was calculated by using perturbation theory, while 
the parameter B represents the contributions from small 
impact parameter and includes both, perturbative and 
non-perturbative parts. The production of bound-free 
electron-positron pairs was calculated also by Bertulani 
and Baur Q who used a semi-classical method in order 
to calculate the production of bound-free pairs at ener- 
gies and for collision systems appropriate for the RHIC 
facility. From these computations it was found that the 
BFPP cross sections for the capture of the electron into 
the ns excited states of the ion decreases with w 1/n 3 , 
which means a factor of 1/8 for the L-shell and to a net 
effect of all ns excited states of approximately 20 % in 
total 

An alternative method was later applied by Rhoades- 
Brown and coworkers Q who performed an " exact" inte- 
gration of the Feynman diagrams by using Monte-Carlo 
techniques. In this work, the cross section for the cap- 
ture of an electron was obtained as the convolution of 
the amplitude of direct and crossed Feynman diagrams 
in Fig. |U 



B(k,q;p ± ) = A<-+\k,q;p x ) 
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with the momentum distribution of the bound-state 
wave function. The results from these Feynman-Montc 
Carlo computations were compared with the Weizsacker- 
Williams calculations by Baur and Bertulani Q and were 
found larger by about a factor of 3, a discrepancy which 
was explained later in a comment by Baur [Io| . In the 
present work, we have calculated the cross section for 
the capture of an electron into the i\-shell by apply- 
ing a Monte-Carlo integration for the lowest-order Feyn- 
man diagrams as shown in Fig . [J This procedure is 
known also in the literature [ll|, [l2| as the two-photon 
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method since the colliding nuclei (nucleus a and nu- 
cleus b) exchange one photon (total two photon) and the 
two-photon-exchange diagrams are proportional to Za. 
In contrast to our previous computations, where plane- 
waves were applied for both the electron and the positron 
fill . fl2j , we here apply bound if -shell wave-functions for 
the electron as well as modified plane-waves functions for 
the positrons that includes a correction due to the dis- 
tortion by the 'bound' electron. In fact, this distortion 
of the positron wave function arises from the necessary 
(re-) normalization of the continuum waves in order to 
account for the reduction of the wave functions of the 
positron near to the nucleus to which the electron is lo- 
calized [HI]. In the literature, these (one-particle) func- 
tions are known as Sommcrfeld-Maue wave functions for 
the positrons and Darwin wave functions for the bound 
electrons. Similar wave functions have been applied also 
in Refs. @, [l3| for studying the captured electrons and 
free positrons. 

In the next section, we first present the formalism for 
evaluating the pair production cross sections with an elec- 
tron bound to one of the ions. Apart from the represen- 
tation of the electron and positron states, this includes 
the analysis and a step- wise simplification of the bound- 
free amplitudes by using the wave functions from above. 
In Section III, then, the differential BFPP cross sections 
are calculated as function of the transverse and longitu- 
dinal momentum, the energy and rapidity, and especially 
for those collision energies of the ions that are relevant 
for the RHIC and LHC facilities. A comparison of our 
Monte Carlo-Fcynman calculations with previous com- 
putations is made. Finally, a few conclusions are drawn 
in Section IV. 



II. THEORETICAL BACKGROUND 

Lowest-order perturbation theory in the framework of 
quantum electrodynamics (QED) has been applied to de- 
rive and calculate the cross section for generating bound- 
free electron-positron pairs in relativistic heavy-ion colli- 
sions. For this electron-positron BFPP process, the sig- 
nature is that the electron is captured by one of the 
colliding ions, while the free positron leaves the colli- 
sion system. In lowest QED order, this process is de- 
scribed by the two Feynman diagrams, the (so-called) 
direct and crossed terms, as depicted schematically in 
Fig. [T] These diagrams represent the leading contribu- 
tions to the bound-free pair production as appropriate 
especially for the high collision energies available at the 
RHIC and LHC facilities. 

For two ions a and b, that collide with high energy, the 
leading contributions to the pair production arise from 
those Feynman diagrams for which each ion interacts ex- 
actly once with the electromagnetic field of the other ion. 
This restriction gives rise to the direct and crossed dia- 
grams from Fig. [TJ and to an gauge-invariant total ampli- 
tude that is Lorcntz covariant. In the sudden (or impulse) 




FIG. 1: Lowest-order Feynman diagrams for the pair pro- 
duction of a bound-free electron-positron pair in heavy-ion 
collisions: (i) direct and (ii) crossed diagrams for the simul- 
taneous capture of the electron into a bound state of target 
(T) ion. In the figure, a and 6 represents the two ions, and q 
is the momentum of the positron Q. 



approximation, the pair creation with simultaneous cap- 
ture of the electron by one of the ions is then described 
by convoluting the electron line in Fig. [T] with the mo- 
mentum wave function of the final bound state In 
the following, we consider the collision of the ions in the 
'collider frame' with zero total momentum of the overall 
system. In these coordinates, the two nonzero compo- 
nents of the vector potential of ion a is given by 



[2nZe] S(q 



A 



d A r ■ 



-sq r 



(4a) 



(4b) 



while, similarly, the potential of ion b is obtained by just 
substituting the impact parameter b — > —b in the expo- 
nent of Eq. (|4a[) and the relative velocity (3 = - — ► — p 
in the formulas above. In these expressions, moreover, 
V refers the velocity of nucleus a moving from right to 
left parallel to the z axis. Nucleus b moves from left to 
right with velocity —v. In this work, we consider only 
the symmetric collisions (equal mass and equal charge) 
of the heavy ions, although the calculations can be done 
easily for the asymmetric collisions. For the overall col- 
lision system, therefore, the classical four-potential 
can be written as [ll|, : 



AH 



(5) 



which describe a retarded Lienard-Wicchert interactions. 

To evaluate the diagrams in Fig. [TJ we need of course 
a proper set of one-particle states in order to represent 
the electron-positron pairs after their generation in the 
field of the moving ions. For the outgoing positron, the 
spinor structure is 



ul+) = 



E, 



2mc 2 



PC 



E { q +) +mc 2 



(G) 
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(for spinors with positive energy Eq + ^ > 0), and where 



Xi/2 denotes a Pauli spinor and s = ±1/2 its spin 



projection. Here, after the creation of electron-positron 
pairs, electron is captured by one of the colliding ions and 
positron becomes free which is described by the plane- 
waves 



*(+) 
g 



AL 



M* u (+) 



(7) 



and together with the (correction) term ^ in order to 
account for the distortion due to the charge of one of the 
nucleus. In expression 0, moreover, 



N+ = e -*n+/2 T (i 



ia+) 



a + 



v+ 



(8) 



is a normalization constant which accounts for the dis- 
tortion of the wave function is acceptable for Za -C 1 
[1, El 03, and where a = e 2 /hc = 1/137 is the fine 
structure constant and v+ the velocity of the positron 
in the rest frame of the ion, into which the electron is 
captured in course of process |T]). For sufficiently large 
energies of the ions, we can approximate v+ = c = 1 by 



using natural units (ft — c — m — 1), and which are used 
throughout this work if not stated otherwise. 

After the pair production (Q]) has occurred, the elec- 
tron is captured by one of the ions and, thus, need to 
be described as a bound state. In a semi-relativistic ap- 
proximation, these electron states are often represented 
by EEl 



= 1 - — a • V non _ rc i (r) , (9) 
Ira 



i.e. in terms of the non-relativistic (ground) state func- 
tion 



1 / Z 



3/2 

*non-rclM = A= ( — ) 



-Zr/a H 



(10) 



of the hydrogen-like ion, and where u represents the 
spinor part of the captured electron and ajj = 1/e 2 the 
Bohr radius of atomic hydrogen. 

Using the positron and electron states from above, the 
direct diagram in Fig. QJi) can be written as: 



(*<->is*i*w) = *EE 

p s 

-EE 



00 dw 
-00 2tt 

dtu 



( + ) 
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(ii) 







/ dh ( 1 H 




/ — 00 V 


2m / 



dV JV+ e-iCp-i)-'' A fc ( r ; M+> - w) 



u I (1 - /3a,) I ng \ (uW | (1 + u ^ 



,(+)\ 



(4 s) - W ) 



(12) 



and where V a and Vj, are the potentials of the two nuclei 
a and b, 



V a = (l-0a,)A° 
V b = (l + f3a z )A° b . 



(13a) 
(13b) 



respectively, and with the vector potentials taken from 
above [cf. Eq. 0]. Apparently, therefore, the overall (di- 
rect) amplitude contains a three-fold integration over the 
coordinates of the bound electron (<i 3 r), the coordinates 
of the free positron (c? 3 r') as well as the integration over 
the frequency u> of the virtually exchanged photons be- 
tween the heavy ions. In the evaluation of Eq. (jll[) . more- 
over, we have used the completeness relation 



E 



y(+) 
Ap 



Y (+) 

Ap 



+ 



Ap 



,(-) 



= 1 (14) 



with Xp being: 



Y (+) 

Ap 



e ipr ui+> 



(15a) 
(15b) 



(s) 

In this notation of the one-particle states, again, Uo-J 

(s) 

refers to the spinor part of the intermediate state of Xp t 
and the summation over the spin and momentum of the 
one particle can be replaced by 



d 3 P 

(2tt) 3 



(16) 



We can perform the integration over r and r' explicitly 
in Eq. (|12p. If we first consider the integral over r, the 
two parts of this integral 
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err ( 1 + —a • V ) * nc 



(r) e ip r 4,(r;w -£<->) 



J d 3 Ta-V^ non . Tel (r)e^ r A a (v;cj-E^) (17) 



can be analyzed independently. Making use of the explicit form of the non-rclativistic Is— function ^non-relC?*) and 
the vector potential A a (r;uj — E^), we can write the first part (on the rhs) of Eq. (fTTj) as 



d 3 rf non _ re i(r)e^ r A (r;a;-^-)) = -[2*Ze]5(po) ^ ***J -= ( — 



W 2 Ze ( Z 



- I — 

3/2 Jipj 



3/2 



d 3 re -Zr/a« 



^" + 7 2 + P ± 



<y(w - - . (18) 



I 

Here, in the second line, we made use of the known inte- Appendix A. In Eq. (|18[) . moreover, po represents the 



gral 



d 3 re ~Zr/a H 



4tt 



(19) 



energy term, p z the longitudinal momentum and, p±(p y ) 
is the transverse momentum of the intermediate state as 
given by Eq. {T] 



(i + p 

together with the Lorentz transformation as displayed in 



The second part of the integral (|17|) 



i 

2m 



d 3 ra- V tf „o„-rciM e*' r ^ a (r;a;- 



3/2 



3/2 



= - J- -L ( — ) 87r 2 Ze a-p*(w- - /?p.) / 



pi 



(20) 



can be evaluated by following similar lines, but it now contains the factor ■ V which arises from the wave function 
of the captured electron. This additional 'derivative' with regard to the coordinates of the electrons can be removed 
by an integration by parts, 



d 3 ra-V e~ Zr/aH 



-in 



= —la. ■ p 



^+P 2 



(21) 



For the overall integral (fTT)) . this gives rise to the expression: 

y°°d 3 r (l + ^a-vj * non _ rel (r)e*- r ^ (r;a;-S(-)) 



a ■ p 



2m 



87r 2 Ze. 



3/2 



6(w - E(~) - 0p z ) l Pj 



(22) 



where E 1 - ' is the energy of the captured electron. 
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Using analogue steps, the integral over r' in Eq. (|12[) can be written as: 



dV N+ e-^P-^-* A b (r ; - lu) 



ipz - q z ) 2 +7 2 (p± - qj.) 2 



-qi)- 



(23) 



where Eq is now the energy of the positron. Thus, by combining both integrals in Eq. (|17| . we obtain for the direct 
BFPP amplitude the explicit expression: 

d 3 p 



iN. 



EE 



(2tt) 3 



dm 
2^ 



Z 
an 



3/2 8(w - EM - 0p g ) r 



7 



pi 



(Pa - Qz) 2 +7 2 (PJ 



q±) 2 
I — 



E, 



(») 
p 



a • p 



2m 



x o7r ze7 



(24) 



(s) 

and where £p is the energy of intermediate state. 

The vector p describes the momentum of the interme- 
diate (electron and positron) states in the field of the ion 
and can be decomposed into its transverse and parallel 
part, p = p± + p z , relative to the motion of the ions. 
In Eq. [2H the two Dirac delta functions gives us the 
component along the velocity of the heavy ions 



E { + ] E(~ 



Pz = 



faz 



2(3 



(25a) 



in terms of the energies and the longitudinal momentum, 
q z , of the positron. Similarly, we can write the frequency 
of virtually exchanged photons as 



w = E^ + f3p z 



E 



(+) 



Pqz 



(25b) 



owing to the conservation of the momentum that is ob- 
tained from the Dirac delta functions in Eq. In 
contrast to the parallel part, however, the transverse mo- 
mentum of the electrons pj_ is not fixed by the kinematics 
of the collision partners, but also depends on the momen- 
tum that is carried by the fields. After integrating the 
Eq. (|24p over lu and p z , and by inserting the values from 
Eqs. ([23)) into Eq. (|24[) . the transition matrix element for 
a fixed spin and momentum state of the positron as well 
as for a given intermediate state can be expressed as 



iN+ 



3/2 



d 2 p± 
(2ny 



r 



1 F(-p x : w„) F(p x - q ± : Lo b ) T q {p ± : +(3) , (26) 



where b is again the impact parameter of the ion-ion 
collision, and the function F(q,w) can be described as 
the scalar part of the field associated with the ions a and 
b in momentum space. Explicit form of these scalar fields 
can be written in terms of the corresponding frequencies 
as 



F(-p ± : L0 a ) 



AirZe 



z- 



Pi 



for the frequency Lu a , and as 



F(P± - q_L : Wb) = 



AnZe^P 2 



K 2 + 7 2/32( p± _ q±) 2 ) 

for the frequency Lu bl respectively 



(27a) 



(27b) 



Owing to the asymmetry in the behavior of the elec- 
trons and positrons in course of the BFPP process, the 
wave functions of the free positron and captured electron 
will differ substantially. This difference gives rise also 
to different expressions for the frequencies of the virtual 
photons as emitted by the two nuclei, i.e. 



E 



(+) 



for ion a, and 



L0 b 



E ( q +) - E(~) - f3q z 



(3 P z 



(28a) 



p(p z -q z ). (28b) 



for ion &, respectively. As mentioned before, these fre- 
quencies are obtained from integrating Eq. (|24|) . Apart 
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from the scalar field of each ion, Eq. (|26j) contains also 
the transition amplitudes T which relates the intermedi- 
ate photon lines to the outgoing electron-positron lines. 
This amplitude depends explicitly on the (relative) ve- 
locity of the ions /?, the transverse momentum p± , and 
the momentum of the positron q, and it is given by 



r q (p± : +/3) 

= EE 



E 



0) 



Ei-)+E 



(+) 



a ■ p 



2m 



x (u |(1 - t 3a z )\ u«) (ug |(1 + 0a z )\ u«) . (29) 

In this amplitude, moreover, the parallel component of 
the intermediate state momentum p z is determined by 
Eq. (|25a|) . Finally, let us note that the integration 



over the impact parameter b in Eq. (|26|) can be car- 
ried out also analytically. Following very similar lines, 
it is possible also to evaluate the cross ed-term amplitude 
(*(-)|S ha |^ +) ) from Fig.ffl 

Having the amplitudes for the direct and crossed dia- 
gram, we are now prepared to write down the cross sec- 
tion for the generation of a free-bound electron-positron 
pair in collisions of two heavy ions 



d 2 b 



q<0 



(-) 



|S|*<+>) 



(30) 



where S = S ab + S ba denotes the sum of the direct and 
crossed terms in Fig. [1] Making use of all the simplifica- 
tions from above, these cross sections for the BFPP can 
be expressed as: 



* = d 2 bJ2 |(* {_) \Sab\ + \S ba \V 



q<0 



/(+) 
■ <? 



4/3 2 



(31) 



with 



and 



A M (q:p ± ) = F(-p ± :u a )F(p ± -q ± :uj b )T q (p ± :+P), 

A [ ~\q : q_L - P_l) = F(p± - q± : uj b ) F(-p± : u) a )T q (q ± - p± : -f3) . 

I 



(32a) 
(32b) 



being some proper products of the transition amplitudes 
and scalar parts of the fields as associated with ions a 
and b. All these functions have been displayed explicitly 
in Eqs. (f2"T|) and (|2"5|) above. Moreover, the square of the 
normalization constant for the positron wave function is 
equal to @, [l3| 



2m, 



(33) 



III. RESULTS AND DISCUSSIONS 

Calculations have been performed for the total pro- 
duction cross sections of bound-free electron-positron 
pairs in relativistic collisions of bare ions. Theoretical 
cross sections are obtained especially for the collisions of 
Au + Au ions at energies relevant for the RHIC facil- 
ity as well as for Pb + Pb ions at LHC energies. These 
cross sections are compared with those for the produc- 
tion of free electron-positron pairs. While, however, the 



Indeed, this constant appears to be very similar to 



2iraZ 



(34) 



TABLE I: Bound- free pair production cross sections ctbfpp 
(in barn) for selected collision systems and cross sections as 
accessible at RHIC and LHC collider facilities. 



i.e. the factor that represents the distortion of the 
positron wave function due to the shielding of the nu- 
cleus by the electron fj~3l. [bH. \vf\ . 



This work Ref. [18] 



RHIC Au + Au at 100 GeV 94.5 94.9 

LHC Pb + Pb at 2957 GeV 202 225 
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-7*100 GeV 
- y=3400 GeV 



10 20 30 40 50 60 70 80 90 



FIG. 2: BFPP cross sections for two different systems as func- 
tions of the nuclear charge Z. BFPP cross sections (in barn) 
for the symmetric collision of bare ions with nuclear charge 
Z at 100 GeV/nucleon (dashed line) and 3400 GeV/nucleon 
(solid line). The capture cross sections increases by about 3 
order of magnitude in going from Z = 10 to Z = 90. 




1500 2000 

T 



FIG. 3: BFPP cross sections for two different systems (Au + 
Au-dashed line and Pb + Pb-solid line) as functions of the 7. 
The magnitude of 7 is going from 10 to 3400. 



free-pair production includes an eight dimensional inte- 
gral, the BFPP cross sections eventually depend only on 
five-dimensional integrals. To evaluate these amplitudes, 
Monte-Carlo techniques were utilized, and the integrands 
have been tested on about 10 million randomly chosen 
'positions' in order to ensure a sufficient convergence of 
our theoretical results. The total numerical errors in the 
computations is estimated to be less or approximately 
five percent. 

As a test of our implementation, first computations 
of the total BFPP cross sections were performed for a 
few selected collision energies. Table U displays these 
BFPP cross sections for the two collision systems from 



above and for those collision energies that arc relevant 
for forthcoming experiments at the RHIC and LHC col- 
lider facilities. Our results for the total cross sections are 
in good-to-excellent agreement with the previous com- 
putations by Meier et al. [l||, especially for the RHIC 
energies of 100 GeV per nucleon. For the much higher 
collision energies of ~ 3000 GeV per nucleon, that will be 
available at the LHC storage ring, our theoretical predic- 
tions for the BFPP cross sections are in contrast lower by 
about 20 %, compared with the computations by Meier 
and coworkers. 

As mentioned above, the correction term was omit- 
ted to the positron wave functions in Eq. {7J, in line 
with previous experience and computations of the free 
(electron-positron) pair production for which a perfect 
agreement with experiment was found by omitting this 
term [ill E3 • We therefore conclude that the distortion 
of the positron states due to the shielding of the electron 
is small and remains negligible for the present computa- 
tions. 

To understand the importance of the bound-free pro- 
cess, Fig. [2] displays the BFPP cross sections for sym- 
metric collisions of ions with charge Z as function of the 
nuclear charge. Cross sections are shown for the two col- 
lision energies E = 100 GeV/nucleon (dashed line) and 
3400 GeV/nucleon (solid line) as important for modern 
accelerators. While the cross sections for the production 
of free pair scale approximately with « (Za)Hn 3 (^i), the 
BFPP cross sections increase with « (Za) s ln(j). In this 
scaling behavior, the reason for an extra factor Z 3 arises 
from the bound wave function of the electron, while an- 
other power in Z comes from the normalization constant 
for the positron wave function. 

Fig-GHshows the total BFPP cross sections for two dif- 
ferent systems as functions of the Lorentz contraction fac- 
tor 7. Our calculation is done in the centcr-of-momcntum 
frame, therefore the relationship between the Lorentz fac- 
tor 7 and the collider energy per nucleon in GeV E/A is 
given by 7 = 1/Vl — V 2 = E/mo, where mo is the mass 
of the nucleons. Results are displayed for Pb + Pb col- 
lisions (solid line) and for Au + Au collisions (dashed 
line). As function of the Lorentz factor 7, the free pair 
production scales with ln 3 (p/) and the bound-free pair 
production with ln(p/). All these results are obtained in 
previous computations [H, [H, EH H3] and, thus, we also 
verify that our calculations gives the similar results. 

Results for the free and bound-free pair production are 
displayed in the figures Q] to [7] within the same graph. At 
collision energies of 100 GeV per nucleon, as relevant to 
the RHIC facility, we have considered Au + Au collisions, 
while Pb + Pb collisions were analyzed at 3400 GeV per 
nucleon, as they are hoped to be reached within the near 
future with the LHC at CERN. The same notation is used 
throughout the figures|3]to[7] Thin-lines and thin-dashed 
lines represent the bound-free pair production and free- 
pair production at RHIC collisions of Au + Au, respec- 
tively. On the other hand, thick lines and thick-dashed 
lines shows the BFPP and free pair production at LHC 
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1Cf 1 10° 10 1 10 z 1Cf 1 10° 10' 10 Z 10 3 

p ± {Transverse Momentum)(Mev/c) P z (Longitudinal Momentum)(MeV/c) 



FIG. 4: The differential cross section as function of the trans- 
verse momentum (p±) of the produced positrons is shown in 
the graph. Calculated differential cross sections are shown for 
the two collision systems Au + Au at 100 GeV per nucleon 
and Pb + Pb at 3400 GeV per nucleon, respectively. When 
compared with the production of free electron-positron pairs, 
the BFPP cross section for the capture of the electron into 
the Is ground state is suppressed by about three orders of 
magnitude for all transverse momenta between 0.1 and 100 
MeV/c. 



FIG. 5: The differential cross section is shown as function of 
the longitudinal momentum (p z ) of the produced positrons. 
The notations are the same as in Fig. [4] 



is observed at LHC energies. For large values of the ra- 
pidity, the bound-free cross sections decay more rapidly 
than for the free production sections. Since the rapidity 
is a function of the energy and momenta, 



collisions of Pb + Pb respectively. 

Fig-LUplots the the differential cross section as function 
of the transverse momentum of the produced positrons. 
From this figure, it becomes clear that the bound-free and 
free pair production distributions display a rather similar 
behavior as function of transverse momentum, although 
the free-pair distribution function is larger by about 3 
orders of magnitude than the BFPP function. Obvi- 
ously, moreover, BFPP distribution function decreases 
much faster with the size of the transverse momentum 
than those for the free-pair production. 

Figs. E]and [5] displays the differential cross sections as 
function of the longitudinal momentum and the energy, 
respectively. Again, cross sections are shown for the free 
and bound-free case. From Fig. Owe find that the ratio 
of the RHIC and LHC increases if the value of the longi- 
tudinal momentum increases. Moreover, since the longi- 
tudinal momentum of the positron is much higher than 
the transverse momentum, the energy of the produced 
positron arises mainly from the longitudinal momentum 
of the positrons. Therefore, the behavior of the differen- 
tial cross section as functions of longitudinal momentum 
and energy must look very similar as seen from Figs. [5] 
and [SI We can co nclude that energy of the positrons, 
E = \pp*\_ + pi + 1, consist of mainly by the longitudinal 
momentum of the positrons. 

Finally, Figj7] plots the differential cross section as 
function of the rapidity. At RHIC energies, the behav- 
ior of the free and bound-free differential cross sections 
of rapidity is almost the same, while a different behavior 



V = 7> ln Po+Pz , (35) 
2 IPO - Pz _ 

the discrepancies in the behavior in the cross sections as 
function of the longitudinal momenta and the energies 
appears closely related to its the behavior as function of 
the rapidity. 




1Cf 1 10° 10' 10 Z 10 3 

P (Energy)(MeV) 



FIG. 6: The differential cross section as function of the energy 
(po) of the produced positrons is shown in the graph. The 
notations are the same as in Fig. U 
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negative value and it is proportional to Z 2 which is ob- 
tained by Bethe-Maximon [22||. So far, these higher-order 
corrections were not included in the computation but we 
plan to derive and calculate this effect in a forthcoming 
work. 
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FIG. 7: The differential cross section is shown as function of 
the rapidity (y). The notations are the same as in Fig. [4] 



APPENDIX A: LORENTZ TRANSFORMATION 
OF THE 4- VECTOR POTENTIAL 



IV. CONCLUDING REMARKS 

In this work, we have investigated the electron-positron 
pair production that is associated with the (simultane- 
ous) capture of the electron into the if— shell of one of 
the ions. Calculations of the cross sections have been per- 
formed especially for two collision systems, Au + Au and 
Pb + Pb, and for energies that are relevant for the RHIC 
and LHC facilities. In the framework of QED pertur- 
bation theory, the lowest-order Feynman diagrams have 
been evaluated by applying Darwin wave functions for 
the Is bound state of the electron and Sommerfeld-Maue 
wave functions for the continuum states of the outgoing 
positron. In line with previous experience, however, we 
have not taken into account the correction term for the 
positron wave function whose influence was estimated to 
be small. We plan to incorporate this term in the future 
to analyze its contribution in further detail. 

Comparison of our theoretical cross sections is made 
with previous computations as far as available. Good 
agreement is found especially for the total cross sections 
of Au + Au collisions at 100 GcV per nucleon as utilized 
at RHIC. For these collision, it is found in particular that 
the free and bound-free pair production cross section be- 
have very similar as a function of energy with an almost 
constant factor of ^ 10 3 with which the bound- free pair 
production is suppressed. This is quite different for colli- 
sion at > 3000 GeV/nucleon, the expected conditions at 
the LHC, where there is discrepancy between BFPP and 
free pair production differential cross sections especially 
for large values of longitudinal momentum, energy and 
rapidity. 

For the total electron-positron pair production, the ef- 
fect of the 'Coulomb correction' is known to play an im- 
portant role 21( due to multi-photon exchange of the pro- 
duced electron-positron with the colliding nuclei. This 
Coulomb correction for the free pair production has a 



In Section II, the integrals (fTSl and [20]) were evalu- 
ated by using the Lorentz-transformcd potentials of the 
Coulomb field of the heavy ions. Here, we derive this 
potential in momentum space and perform the Lorentz 
transformation onto it. In the rest frame of an ion with 
the nuclear charge Z, fixed to the coordinates (0,b/2,0), 
the four-vector potential is given by: 



A' 







A Q = 



-Ze 



c' 2 + (y' -6/2) 2 + z' 2 ]V2 



(Al) 



In momentum space, this vector potential can be ex- 
pressed by 

A' (q) = J d^xe^'^A^x) 

= -[2irZe]S(q )e [ - iq y% ] 



d 6 r 



and with 



x +(y -b /2) + z 



(A2) 



(A3) 



Moreover, since A = and A^ must transforms like a 
four-vector, we obtain 



Ml) 


= 7(A) - 


(3A[) 


= 7A) 


Mo) 


= 7(4 - 




= -7/?A) 




= a', = 


0. 





(A4) 

Thus, by doing the integration over r' on the rhs of 
Eq. (|A2[) and by performing the Lorentz transformations, 
we obtain the potential for the colliding ions: 



A = - [8^ 2 Ze](5(<7o + /3gi)7 2 



(A5) 
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and where q\ refers to the longitudinal momentum, 
the transverse momentum and go is the energy term. 
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